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ABSTRACT 

We have used Virtual Observatory technology to analyse the disk scale length and central surface 
brightness hq for a sample of 29955 bright disk galaxies from the Sloan Digital Sky Survey. We use the 
results in the r-band and revisit the relation between these parameters and the galaxy morphology, 
and find the average value (/lo) = 20.2 ±0.7 mag arcsec -2 . We confirm that late type spirals populate 
the lower left corner of the r^—jiQ plane and that the early and intermediate spirals are mixed in this 
diagram, with disky ellipticals at the top left corner. We further investigate the Freeman Law and 



confirm that it indeed defines an upper limit for /io in bright disk galaxies with 



mag 



< 17.0, and that 



disks in late type spirals (T > 6) have fainter central surface brightness. Our results are based on a 
volume corrected sample of galaxies in the local universe (z < 0.3) that is two orders of magnitudes 
larger than any sample previously studied, and deliver statistically significant implications that provide 
a comprehensive test bed for future theoretical studies and numerical simulations of galaxy formation 
and evolution. 

Subject headings: Galaxies: formation; Galaxies: structure; Galaxies: evolution 



1. INTRODUCTION 

The current mainstream galaxy formation paradigm 
states that galaxy disks form within dark matter halos 
and that there is an intimate relation between the scale 
length rd of the disk and that of the halo. The mass 
distribution of the disk is entirely set by the r<j and for 
example in the exponential case, 60% of its total mass 
is confined within two scale lengths and 90% within four 
scale lengths. Moreover, the angular momentum of the 
disk is set by r^ and the mass distribution of its host 
halo, and the fact that the angular momentum vectors 
are aligned suggests that there is a physical relation be- 
tween the two. During the formation process, galaxy 
mergers and associated star formation and feedback pro- 
cesses play a crucial role in the resulting structure, how- 
ever, the observed sizes of disks suggest that the com- 
bination of these physical processes yield that galactic 
disks have not lost much of the original angula r momen- 
tum acquired from cosmological torqu es (e.g., White & 
Rees ||1978| |Fall fc Efstathiou ||l980l The hierarchical 
and iritall models predict com parable r^ and in a cold 
collapse scenario ( |Vand [1959 ) since angular momentum 
is conserved, immediately after the collapse the gas is 
supported by rotation so that it quickly collects in the 
region where a disk forms with substantially higher ro- 
tation velocity than the halo is to form. A large r^ disk 
forms when the disk mass is smaller than the halo mass 
over the disk region, and vice versa, a small disk forms 
when the mass of the disk dominates the mass of the halo 
in any part of the disk. The self gravitating disk will also 
modify th e shape of the rotation curve nea r the centre of 
a galaxy (Gelato & Sommer-Larsen 1999) and the disk 
is then set to undergo secular evolution. The natural im- 
plication of this scenario is that the r^ dictates the life of 
a disk, and consequently, is a prime factor which deter- 
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mines the position of a galaxy on the Hubble sequence. 

One prominent indicator for a smooth transition from 
spiral toward SO and disky ellipticals is provided by the 
rd-fio diagram where /io is the central surface brightness 
of the disk, where spirals and SOs are mixed and disky 
ellipticals populate the upper left cor ner of this diagram 
( |Kent |[l985| |Scorza fc Bender p~995|). Anot her instruc- 
tive relation is the Freeman Law ( h reeman 1970) which 
relates (j,q to the galaxy morphological type. Although, 
some studies have found that the Freeman Law is an arte- 
fact due to selection effects 



flT98lj IScorzalT van den Bosch ||1998p , recent work have 




et al. 



shown that proper consideration of selection effects can 



be combined with kinematic studies to explore an evolu- 
tionar y (e.g .. van der Kruit~||1987 de Jong [1996 Koda 
2000| jvan der Krmt*P 002 ) 



In the comparison between theory and observations, 
two issues complicate matters. On the theory side map- 
ping between initial halo angular momentum and is 
not trivial, partly due to the fact that commonly the ini- 
tial specific angular momentum distribution of the vis- 
ible and dark component favour disks which are more 



centrally concentrated disks than expone ntial (e.g., de 



Jong fc Lacey |2000 van den Bosch |200T ) . Observation- 



ally, comprehensive samples have yet not been studied 
and the mixture of different species such as low and high 
surface brightness galaxies complicate the measurements 
of disk para meters ( McGaugh et al. |1995 Graham & de 



mok1|200T| ) 



Here we analyse the ra and fio from an unprecedent- 
edly large sample of bright disk galaxies in the nearby 
universe (z < 0.3) using the Sloan Digital Sky Survey 
(SPSS) Data Re lease 6 (|York et al. ||2000| |Adelman- 
i McCa rthy et al. 
(|20l0| (hereaft" 



F 



SPSS) Data Release 6 ( |York et al 
. As detail 
< r ±'10), both parameters were robustly 



2008). As detailed in IFathi et a. 



determined for 30374 galaxies in the r-band (only 29955 
used here as described in section [2| , whereas in other 
SPSS bands, the derived parameters are valid only for 
subsets of this sample. In the g, i, and z-band (« 27000- 
30000 galaxies) the sample sizes are comparable to the 
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Fig. 1. — Comparison of our r-band and ^iq with |Gadotti | 
( |2009| > , where the ratios between the two studies are shown in the 
corresponding panel (y-axis). 

one presented here. In the w-band, and /io were ro- 
bustly derived for a few hundred objects and therefore, 
not considered here. Throughout this paper, unless oth- 
erwise stated, we use disk parameters in the r-band and 
investigate the two relations mentioned above, in order 
to provide a comprehensive test bed for forthcoming cos- 
mological simulations (or analytic/semi-analytic models) 
of galaxy formation and evolution. 

2. THE DERIVED PARAMETERS AND SAMPLE 
CORRECTION 

The derivation of ra and /io for the SDSS sample is de- 
tailed in F10, and testing the code with standard IRAF 
routines delivered fully consistent results. Comparing 
disk parameters with those from several published sam- 
ples jde Jong |1996[[Gr"aham |2003[ |Anderson et al.|2004 



Koopmann et al.||2006| |Balcells et al.||2007p failed to de 



liver any objects in common. However, cross matching 
with the multi-band image fitting and multi-component 
(disk, bar, bulge, and central source) profile d ecompo- 
sition of nearly 1000 galaxies based on SDSS ( Gadotti 
2009) returned 322 disk galaxies for which we derived 
comparable values (see Fig. [T]). As outlined in F10, we 
made use of the Sky View services to obtain cut-outs from 
the SDSS tiles, with the galaxy counts derived in the 
reference frame of the tile that is at the centre of the 
Sky View cut-out. The disk central surface brightness 
can thus be derived using the procedure provided by the 
SDSS services. Applying pixel area correction and Tol- 
man cosmological dimming correction (z + l) 4 we were 
able to obtain the //o in correct units of mag arcsec -2 . 
The longitudes provided by the SDSS services show that 
for our sample ||6|| > 20 °, and we use the stand ard ex- 
tinction correction law of Giovanelli et al. ( 1994[) to cor- 
r ect for this effect. Following the description in | valentijn| 
(1990) and de Jong | ( |1996[ ), we use a range of internal 
extinction parameters (0 for optically thick galaxies to 1 
for transparent galaxies) for every diagram to quantify 
this effect throughout our analysis, but present diagrams 
calcul ated with the extinction parameter of 0.26 (de Jong 
p96f. 

We also derived the asymmetry parameter and the cal- 
culation of the inverse concentration parameter (F10). 
By examining a small subset of galaxies for which reli- 
able morphologic al classifications are available form the 
LEDA catalogue (Paturel et al. 2003 ), we assert that the 



n 



stellar velocity dispersion, asymmetry, and concentration 
can be used with caution as indicators of the position of 
a galaxy along the Hubble sequence. We stress, how- 
ever, the data points have large spread and at best have 
a correlation coefficient of determination 9\ z = 0.31, well 
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Fig. 2. — The r^-fio plane of our sample (density plot of the 
29955 galaxies), with the 282 reliably classified galaxies illustrated 
as coloured points where the colour indicates the numeric morpho- 
logical type T. The dashed line (slope 2.5) illustrates the line of 
constant disk luminosity, and the crosses at the top rig ht illustrate 
the seven disky ellipticals fromlScorza & Bender I d 1 995 b - At the top 
and right, the normalised distribution ot each parameter is shown 
before (dotted) and after (solid) applying the volume correction 
described in section [2] 



below one-sigma level (see F10 for details). 

Our derived and /*o for the SDSS sample will have to 
be corrected for the frequency of galaxies with a certain 
property in a volume to be included in the sample. The 
volume over which we can detect low luminosity galax- 
ies is smaller than the volume over which we can de- 
tect high luminosity galaxies. In order to calculate the 
number density contribution of each galaxy, we calculate 
the maximum volume over which it is observable Vmax- 
We follow the volume weighting recipe described in (e.g. 



van der Kruit 



1996), and weight each 



where Vtnt is the 



1987 |de Jong 
galaxy by the calculated l/ to t/ 
volume contained w ithin the distance to farthest galaxy 
in the sample (e.g., Davies |[l990| |Blanton et al. |2 003). 



To calculate the volumes, we assume a ACDM Universe 
with the Hubble parameter Hq = 70 km s -1 Mpc -1 , 
SIa = 0.7, and Qm = 0.3, and for this calculation use 
g-band luminosities to calculate luminosity distances as 
these are more representative of B-band values. Chang- 
ing the cosmology has no significant effect on the analysis 
carried out here. 

This method has the disadvantage that if the near part 
of the sample is overdense, where a magnitude-limited 
survey is most sensitive to low-luminosity galaxies, then 
the Vtot/Vmax w iU become very large for the faintest ob- 
jects. To avoid this problem, we make a cut in the sam- 
ple by removing galaxies with unreasonably large weights 
(> 10000). This leaves us with 29955 galaxies at z < 0.3 
and r mag < 17.0, for which we obtain volume corrected 
parameter distributions (see e.g., Fig. [2|. 



Revisiting the Scale Length-/io Plane and the Freeman Law 
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Fig. 4. — The r^—no plane colour coded by asymmetry (top), 
concentration parameter (middle), and velocity dispersion a (bot- 
tom), where the ranges have been chosen to amplify the effects 
discussed here. 



3. THE r d -n PLANE 

Figure [2] shows that the derived and /x follow 
the linear trend previously discussed i n the context of 
the formation and evolution of disks (de Jong 1996 



Courteau |1996{|van den Bosch |1998"l|Graham &: de Blok| 



|2001 ). We find that fj,o has a Gaussian distribution with 
(Ho) = 20.1 ± 0.7 mag arcsec -2 with a linear trend seen 



in Fig. [5] (applying different internal extinction parame- 
ters changes this mean value by 0.2 mag arcsec -2 ). The 
top right corner is enclosed by the constant disk lumi- 
nosity line, void of objects. The top right corner is also 
the region where the disk luminosity exceeds 3L*, thus 
the absence of galaxies in this region cannot be a selec- 
tion effect since big bright galaxies cannot be missed in 
our diameter selected sample. However, it is clear that 
selection effect plays a role in populating the lower left 
corner of this diagram. Analogue to the r^-fiQ plane, the 
Tully-Fisher relation implies lines of constant maximum 
speed a disk can reach. 
The type d ependency of the r^-zzp plane was previously 



presented by Graham fc de Blok 
high and low surface brightness ga" 



(2001) who used 120 



type-dependency of the diagram. 



axies to explore the 
Our 282 well-classified 



galaxies (illus trated with coloured dots in Fig. 21) follow 
the results of |Graham & de Blok I (120011) and Gadotti 



( 2009 1 . Disks of intermediate and early type spirals have 



higher fiQ while the late type spirals have lower /xo, and 
they populate the lower left corner of the diagram. The 
r,i, on the other hand, does not vary as a function of 
morphological type (F10). 

Figure [3] shows that a forth quantity is equally impor- 
tant in this analysis. T h e tota l gal axy mass obtained 
from [Kauffmann et al. I ((2003 ) and Brinchmann et al. 



[] ( 2004[ ), separates the data_along lines parallel to the 
dashed li 



lines drawn in Fig. [2] This is indeed also con- 
firmed by the Tully-Fisher relation. Moreover, we val- 
idate that the lower mass galaxies are those with type 
> 6. 

Our data also show the expected trends: Asymmetry 
and concentration increase for later types, and central 
stellar velocity dispersion decrease for later type spiral 
galaxies, however, we note that these correlations are 
well below one-sigma confidence level (see end of sec- 
tion [2] and F10). In Fig. H we illustrate r^-fJ-o plane 
where the data points have Deen colour coded by asym- 
metry, concentration, and velocity dispersion as a rough 
indication of the galaxies in different types. 

The top panel of Fig. ^indicates that the higher asym- 
metry galaxies populate a region more extended toward 
the bottom right corner with respect to the high asymme- 
try galaxies. The middle panel shows an opposite trend, 
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Fig. 5. — Volume corrected distribution of disk central surface 
brightness fio versus morphological type T, asymmetry, concen- 
tration parameters, and stellar velocity dispersion. The red dots 
illustrate the 282 galaxies for which robust morphological types 
were given in the LEDA catalogue (type error smaller than 0.5) 
and the open circles show the average fio f° r each type. 



and the bottom panel shows that larger velocity disper 
sion has the same effect as a symm etry. In full agreement 



We plot the Freeman (1970) diagram, where the den- 
sity plot indicates that disks have constant surface 
brightness put to Scd galaxies ( types T = 5, top panel 
of Fig. [5]). The volume corrected density plot also shows 
that the vast majority of the sample is of morphological 
type earlier than 5, thus introducing a bias in the ob- 
served trend. When we consider only the galaxies with 
robust morphological classification (red dots in the top 
panel of Fig.j5l), we n ote that /x n dro ps for galaxies of 
T > 6 (c.f., Fig. 3 of |de Jong |[1996| > with the signifi- 
cant advantage that our sample contains five times more 

f alaxie s than the number of late type spirals of |de Jon"g| 
(1996). Furthermore, using the entire sample of 29955 
galaxies, we confirm that the Freeman Law indeed defines 
an upper limit for /iq in bright disk galaxies at z < 0.3 
(low-surface brightness galaxies are not present in our 
sample). Combined with our previous results, i.e. that 
Td varies by two orders of magnitude independent of mor- 
phological type, this result implies that disks with large 
scale lengths (and typically higher mass) not necessarily 
have higher /xq. 

We use the asymmetry, concentration, and stellar ve- 
locity dispersion as indicators of morphological type. 
The lower three panels of Fig. [5] illustrate that the den- 
sity plots for the entire sample and the galaxies. None of 
the trends are statistically significant in exploring these 
parameters against the Freeman diagram. The lack of 
strong correlation may well be attributed to the fact that 
three parameters that we have used here do not corre- 
late strongly with morphological type. The concentra- 
tion para meter is further complicated by the nature of 
the bulges Gadotti ( 2009 ) . We conclude that these pa- 
rameters cannot be used to further explore the Freeman 
Law. 

5. CONCLUSIONS 

We have derived and analysed r^ and /^o in the r-band 
from SDSS images of the largest ever volume corrected 
sample of disk galaxies in the local universe (z < 0.3 
and r mag < 17.0). We confirm that early and inter- 
mediate type spirals are mixed in the r^-^o plane and 
late type spirals (T > 6) populate the lower left cor- 
ner of this plane. Varying the internal extinction value 
does not change the trends found here, and the average 
value (fio) = 20.1 ±0.7 mag arcsec -2 only varies slightly. 
Moreover, the Freeman Law defines an upper limit for 
jtzrj, where disk galaxies of type T > 6 have fainter /iq. 

These results in r-band are comparable with those 
in other SDSS bands, and the average no in other 
bands are: g-band (20.3 ± 0.7 mag arcsec -2 ), i-band 
(20.1±0.7 mag arcsec -2 and full y consistent with the me- 
dian value from |GadottT| ( |2009[ ) ) , and z-band (21.6±0.7 
mag arcsec -2 ). Given the intrinsic scatter in each band, 
we do not find these differences significant to inter- 
pret the effects of du st and/or stellar populations (c.f., 



iPeletier et al. 19941 moreover that presented in FIO. 



with |Graham fc de Blok~| ( |2001[ ), Fig. [| confirms that Furthermore, our results imply that disks with large scale 

lengths (and typically higher mass) not necessarily have 

brighter fiQ. 

As earlier discussed by e. g., Kent (|1985[); van der Kruit 
H( |l987| ); |Scorza fc Bender | ( p9l^|van der Kruit | fe0(ffl 
and many others, selection effects could cause artificially 
tight or incorrect correlations. In the two relations anal- 
ysed here, we find typically larger scatter than previ- 



galaxies with early and intermediate morphologies are 
mixed along a linear slope of 2.5 in the r^-^o plane, co 
inciding with the region populated by SOs as shown by 



Kent 



Bender 



( 1985 1 and disky ellipticals shown by Scorza & 



4. THE FREEMAN LAW 
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cms analyses, and although our sample represents bright 
disks, the sample size adds credibility to our findings. 
These results are fully consistent with the common un- 
derstanding of the rd-^o plane and the Freeman Law, 
and they contribute to past results since they are based 
on a sample which is two orders of magnitude greater 
than any previous study, with more than five times more 
late type spiral galaxies than any previous analysis. 
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